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ABSTRACT 

We study the accuracy with which weak lensing measurements could be made from a future 
space-based survey, predicting the subsequent precisions of 3-dimensional dark matter maps, 
projected 2-dimensional dark matter maps, and mass-selected cluster catalogues. As a baseline, 
we use the instrumental specifications of the Supernova/ Acceleration Probe (SNAP) satellite. We 
first compute its sensitivity to weak lensing shear as a function of survey depth. Our predictions 
are based on detailed image simulations created using 'shapelets', a complete and orthogonal 
parameterization of galaxy morphologies. We incorporate a realistic redshift distribution of 
source galaxies, and calculate the average precision of photometric redshift recovery using the 
SNAP filter set to be Az — 0.034. The high density of background galaxies resolved in a wide 
space-based survey allows projected dark matter maps with a rms sensitivity of 3% shear in 1 
arcmin^ cells. This will be further improved using a proposed deep space-based survey, which 
will be able to detect isolated clusters using a 3D lensing inversion techniques with a la mass 
sensitivity of approximately 10^"^ Mq at z = 0.25. Weak lensing measurements from space will 
thus be able to capture non-Gaussian features arising from gravitational instability and map out 
dark matter in the universe with unprecedented resolution. 

Subject headings: dark matter — gravitational lensing — large-scale structure of universe — space 
vehicles. 
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1. Introduction 

Weak gravitational lensing lias now been estab- 
lished as a powerful technique to directly mea- 
sure the large-scale mass distribution in the uni- 
verse (for reviews, sec Mellier 1999; Bartelmann & 
Schneider 2001; Refrcgicr 2003). Several groups 
have measured the coherent distortion of back- 
ground galaxy shapes around known galaxy clus- 
ters {e.g. Joffre et al. 2000; Dahle et al. 2002) 
and also statistically in the field {e.g. van Waer- 
beke ct al. 2001; Bacon et al. 2002; Hoekstra et 
al. 2002; Jarvis et al. 2003). Ever-growing sur- 
veys using ground-based telescopes are beginning 
to yield useful constraints on cosmological param- 
eters (Bacon et al. 2003; Brown et al. 2002; Hoek- 
stra, Yee, & Gladders 2002; van Waerbeke et al. 

2002) . The first two clusters selected purely by 
weak lensing mass have now been found and spec- 
troscopically confirmed by Wittman et al. (2001, 

2003) . 

Weak lensing is of such great interest for cos- 
mology because it is directly sensitive to mass. 
Other observations have traditionally been lim- 
ited to measuring the distribution of light and 
linked to theory via complications like the mass- 
temperature relation for x-ray selected clusters 
(PierpaoU, Scott & White 2001; Viana, NichoU & 
Liddle 2002; Huterer & White 2003) or the ubiqui- 
tous problem of bias (Weinberg et al. 2003). Weak 
lensing measurements first avoid these problems, 
then have even been used to calibrate other tech- 
niques (Huterer & White 2003; Gray et al. 2002; 
Hoekstra et al. 2002b: Smith et al. 2003). The 
high resolution, galaxy number density and sta- 
ble image quality available from space-based weak 
lensing data will allow maps of the projected dis- 
tribution of dark matter to be reconstructed at 
unprecedented resolution. The mass power spec- 
trum can be sliced into multiple redshift bins us- 
ing photometric redshifts: providing a long lever 
arm for constraints on the evolution of cosmo- 
logical parameters. Even three-dimensional mass 
maps, marginally feasible from the ground (Bacon 
& Taylor 2002), are likely to be sensitive to over- 
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densities as small as galaxy groups from space. 

Mass-selected cluster catalogs can also be ex- 
tracted from such maps (Weinberg & Kamionkowski 
2002, Hoekstra 2002). Cluster counts, and the 
quantitative study of high-sigma density pertur- 
bations or higher order shear correlation func- 
tions (Bernardeau, van Waerbeke & Mellier 1997; 
Cooray, Hu & Miralda-Escude 2000; Munshi & 
Jain 2001; Schneider & Lombardi 2003) are one of 
the most promising routes to breaking degenera- 
cies in the estimation of cosmological parameters 
including and w, the dark energy equation 
of state parameter. Furthermore, studying well- 
resolved groups and clusters individually, rather 
than statistically, will lead to a better understand- 
ing of astrophysical phenomena, the nature of 
dark matter and the growth of structure under 
the gravitational instability paradigm {e.g. Dahle 
et al. 2003). 

In this paper, we predict the general sensitiv- 
ity to weak lensing of a space-based wide field 
imaging telescope, taking as a baseline the spec- 
ifications of the proposed Supernova/ Acceleration 
Probe (SNAP) satellite. Instrument characteris- 
tics, including the PSF, ellipticity patterns and 
image stability have been studied by Rhodes et 
al. (2003; paper I). In §2 we introduce detailed 
simulated images that have been developed us- 
ing shapelets, an orthogonal parameterization of 
galaxy shapes (Massey ct al. 2003, Refregier 2003). 
The simulated images contain realistic populations 
and morphologies of galaxies as will be seen from 
space, modelled from those in the Hubble Deep 
Fields (HDFs; Wilhams et al. 1996, 1998). These 
shapclct-galaxies can be artificially sheared to sim- 
ulate gravitational lensing. The subscriuent re- 
covery accuracy of the known input shear is dis- 
cussed in §3. We discuss the accuracy of SNAP 
photometric redshifts in §4. Those two measure- 
ments are combined to predict the accuracy of pro- 
jected dark matter maps, 3-dimensional dark mat- 
ter maps and mass-selected cluster catalogues in 
§5. We draw conclusions in §6. Our results are 
used the predict the accuracy of cosmological pa- 
rameter constraints in Refregier et al. (2003; paper 
HI). 
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2. Image Simulations 

In this section, we describe our method for sim- 
ulating realistic images, aiming to closely resem- 
ble images observed with a space-based telescope. 
These simulations are part of a full pipeline which 
allows us to propagate the effects of perturbations 
in the instrument design onto shear statistics and 
cosmological paramcitcrs. Example simulated im- 
ages are shown in figure 1. 

2.1. Procedure 

The shapelet formalism (Refregier 2003; Re- 
fregier & Bacon 2003; summarized in §2.2) has 
been used to model all the galaxies in the HDFs. 
Using just a few numbers, this parameterization 
captures the detailed morphology of the galaxies, 
including spiral arms, arbitrary radial profiles and 
irregular substructure. The parameters for each 
galaxy are stored in a multidimensional parameter 
space. This is then randomly re-sampled, to simu- 
late new and unique galaxies with realistic proper- 
ties as compared to those in the original HDFs. A 
detailed description of the simulation procedure 
and performance can be found in Massey et al. 
(2003). 

The simulated images are built up with galax- 
ies of all types (spiral, elliptical and irregular) in 
their observed proportions, with realistic number 
counts and a size distribution reproducing that 
in the HDF. Their morphology distribution as a 
function of magnitude also reproduces that in the 
HDF. Most importantly, all of these objects pos- 
sess a precisely known shape, magnitude, size and 
shear. The amount of shear can be adjusted in 
shapelet space as an input parameter. 

Observational effects including PSF convolu- 
tion, pixelization, noise and detector throughput 
are then incorporated in the simulations. In §2.3 
we describe the engineering specifications we have 
used to emulate the performance of the SNAP 
satellite. In §3 we then attempt to recover the 
known input shear from these realistic, noisy im- 
ages using existing (and independent) shear mea- 
surement methods. 

2.2. Shapelets 

Here we briefly describe the idea of shapelets, 
which is at the core of our image simulation pack- 



age. More comprehensive details are available 
in Refregier (2003), Refregier & Bacon (2003), 
and Massey et al. (2003). Shapelets are an or- 
thonormal basis set of 2D Gauss-Hermite fimc- 
tions. They can be used to model any localized ob- 
ject by building up its image as a series of succes- 
sive basis functions, each weighted by a "shapelet 
coefficient" , rather like a Fourier or wavelet trans- 
form. Each polar basis state and shapelet coef- 
ficient can be identified by two integers: n > 
describing the number of radial oscillations, and 
m e {—n,n} the azimuthal oscillations, or rota- 
tional degrees of symmetry. The basis is complete 
when the series is summed to infinity, but it is 
truncated in practice at a finite Umax- This of- 
fers image compression because an object is typi- 
cally well-modelled using only a few shapelet co- 
efficients. 

Conveniently, the shapelet coefficients are 
Gaussian-weighted multipole moments (with the 
rms width of the Gaussian known as the shapelet 
scale size /?), as commonly used in various as- 
tronomical applications. The n = 2 states are 
thus Gaussian-weighted quadrupole moments, the 
n = 4 states octopole moments, etc. Shapelet ba- 
sis functions also happen to be eigenstates of the 
2D Quantum Harmonic Oscillator, with n and m 
corresponding respectively to energy and angular 
momentum quantum numbers. This analogy sug- 
gests a well-developed formalism. For instance, 
shears and dilations can be represented analyti- 
cally as a or ladder operators (Refregier 2003); 
and PSF convolutions as a trivial bra-ket matrix 
operation (Refregier & Bacon 2003). 

Massey et al. (2003) demonstrate how HDF 
galaxies can be represented as shapelets and 
then transformed by slight adjustments of their 
shapelet coefficients into new shapes. This pro- 
cess produces genuinely new but realistic galaxies, 
as proved by the similar distributions in HDF 
and simulated data of commonly used diagnos- 
tics from SExTRACTOR (Bertin & Arnouts 1996) 
and galaxy morphology packages {e.g. Conselice, 
Bershady & Jangren 2000). 

2.3. SNAP Simulations 

For this work, our image simulations have 
been tuned to the instrument and specifications 
of the proposed SNAP mission (paper I; Alder- 
ing et al. 2002; Kim et al. 2002; Lampton et al. 
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Fig. 1. — 30" X 30" portions of simulated SNAP 7-band images, using tlie PSF sliown in figure 2. Left panel: to the 
deptli of tiie proposed SNAP wide survey. Right panel: to the depth of the Hubble Deep Fields. The SNAP deep 
survey will be some 2 magnitudes deeper than the latter, but further real data from the ACS on HST are needed to 
simulate this depth accurately. 



2002a, 2002b; Pcrlmuttcr ct al. 2002). The SNAP 
strategy includes a wide, 300 square degree sur- 
vey (with 4 X 500s exposures reaching a depth of 
AB 27.7 in R for a point source at 5(t), and a 
deep, 15 square degree survey (120x4x300s to 
AB 30.2). For an exponential disc galaxy with 
FWHM=0.12", these limits become 26.6 and 28.9 
respectively. 

The predicted SNAP PSF at the middle of the 
illuminated region of the focal plane is illustrated 
in figure 2. Following the analysis of paper I, this 
was obtained for the current satellite design, us- 
ing raytracing, aperture diffraction and CCD dif- 
fusion. In this paper we also illustrate the de- 
composition of the SNAP PSF into shapelets. As 
shown on the top panel of figure 2, our model in- 
cludes the second diffraction ring and is accurate 
to nearly one part in 10^. It does not include much 
of the extended low-level diffraction spikes, which 
we ignore. Convolution with this residual PSF 
pattern adds less than 0.7% to the ellipticity of 
any exponential disc galaxy that passes the size 
cut into the lensing catalog (see §3.2). Given the 
further factor of in equation [4], to convert 
ellipticity into shear, this residual thus has a neg- 
ligible impact upon shear measurement within the 



accuracy of the current methods. 

Simulated images used to calibrate the shear 
measurement method (see §3.2) were first sheared 
and then convolved with the full SNAP PSF shown 
in figure 2. For this application, it is essential 
that the shearing is applied before the smearing, 
just as occurs in the real universe. Shear mea- 
surement methods have been designed to correct 
for precisely this sequence of events. However, our 
simulated galaxies were modelled on real HDF ob- 
jects which had already been naturally convolved 
with the WFPC2 PSF when the HDF images were 
taken. Consequently, our simulated objects in 
§3.2 exhibit smoothing from both a circularised 
WFPC2 PSF, (plus shearing), plus a SNAP PSF. 
This double PSF artificially reduces the rms ellip- 
ticity of galaxies by approximately ~2% and in- 
creases the size of a point source by 22%. One 
should note that the first PSF convolution occurs, 
and the galaxy orientations are randomized, all 
before shearing. This effect therefore corresponds 
to a small alteration in the intrinsic shape distri- 
bution of galaxies but does not bias the shear mea- 
surement (see discussion in Massey et al. 2003). 

Simulated images used to predict the lensing ef- 
ficiency as a function of exposure time (see §3.3) 
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Fig. 2. — Shapelet decomposition of the proposed 
SNAP PSF. Top panel: a horizontal slice through 
the center of the real (solid line) and shapelot- 
reconstructed (dashed line) PSF. The middle panels 
shows, in 2-dimensions, the real PSF, its recovery us- 
ing shapolets and the residual difference between the 
two, from left to right respectively. Bottom panel: the 
moduli of the corresponding polar shapelet coefficients 
with order up to Umax = 12. Note that all intensity 
scales are logarithmic. The circular (m = 0) core is 
modelled to an accuracy of about 10~'' and the begin- 
nings of six-fold symmetric structure is seen as power 
in the m = ±6, ±12 shapelet coefficients. 



were produced differently. For these, we needed 
to ensure realistic size distributions and number 
counts in the simulations. The galaxies had no 
artificial shear added: they just have a scatter of 
ellipticities due to their own intrinsic shapes. We 
convolved these galaxies by the PSF difference be- 
tween the HST and SNAP. This is obtained by de- 
convolving the WFPC2 PSF from the SNAP PSF 
model, in shapelet space. Smoothing an object 
with this smaller kernel is enough to convert it 
from an observation with HST to one with SNAP, 
although without inputting shear. 

Example simulated images are shown in figure 
1 for the wide SNAP survey (left panel) and to 
the depth of the HDF (right panel) . They include 
a noise model consisting of both photon counting 
error and a Gaussian background. These compare 
well with real deep HST images (see Massey et al. 
2002). The SNAP deep fields wiU be about 2 mag- 
nitudes deeper than the HDF. However, deeper 
surveys with the ACS on board HST are awaited 
to accurately model galaxies at this depth. Figure 
3 shows the size-magnitude distribution of the sim- 
ulated images to both depths (top panels). Again, 
the simulations reproduce the statistics of the real 
HDFs (bottom panels). 

2.4. Limitations of the Simulations 

The SNAP wide survey strategy includes four 
dithered exposures at each pointing. This will 
enable the removal of cosmic rays and, if neces- 
sary, the simultaneous measurement of instrumen- 
tal distortions. Because of the high orbit and slow 
thermal cycle, instrument flexure and the PSF are 
expected to be very stable (sec paper I). It should 
therefore be possible to map internal distortions 
and compensate for them even on small scales, us- 
ing periodic observations of stellar fields. Conse- 
quently, neither cosmic rays nor astrometric dis- 
tortions are added to the simulations. 

The SNAP CCD pixels are 0.1" in size and thus 
under-sample the PSF. To compensate for this, 
the dithered exposures will be stacked, as usual 
for HST images, using the DRIZZLE algorithm 
(Fruchter & Hook 2002). Alternatively, galaxy 
shapes may be fitted simultaneously from several 
exposures. DRIZZLE recovers some resolution, 
and will be particularly effective for the multiply- 
imaged SNAP deep survey, but has the side-effect 
of aliasing the image and correlating the noise in 
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adjacent pixels. Wc have not yet included this 
entire pipeline in the simulations, but merely im- 
plemented a smaller pixel scale and model back- 
ground noise that is higher in each pixel (although 
uncorrelated). Following the example of the Hub- 
ble Deep Field final data reduction, we choose 
0.04" pixels. Unfortunately, the detection and 
shape measurement of very faint galaxies is sen- 
sitive to the precise noise properties of an im- 
age. Because of these instabilities, our simulated 
images are only reliable down to approximately 
I ~ 29.5 (see Massey et al. 2003). This is just be- 
low the magnitude cut applied by our shear mea- 
surement method at / = 29.1. A further investi- 
gation will include full use of DRIZZLE and more 
detailed noise models. This will also address the 
issue of pointing accuracy, and consider the con- 
sequences of 'dead zones' around the edges of the 
pixels which house the CCD electronics and are 
therefore unresponsive to light. 

SNAP wide simulated image SNAP deep(er) simulated image 
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Fig. 3. — Size vs magnitude as determined by SEx- 
TRACTOR with a S/N cutoff at u = 1.5. Top panels 
are for simulated SNAP / band images of the same 
size as the Hubble Deep Field. For reference, the bot- 
tom panels are of the HDFs themselves using the same 
SEXTRACTOR parameters. 

The image simulations are based upon the 
galaxies in the HDF, which is itself a special re- 



gion of space selected to contain no large or bright 
objects. As a result, our simulations do not yet 
include these either. The source catalog is being 
expanded as GOODS ACS data becomes publicly 
available. 

The image simulations are currently mono- 
chromatic, in the HST FSIAW (hereafter /) filter. 
Since gravitational lensing is achromatic, shear 
measurement can be performed in any band: in- 
deed, all tested shear measurement methods so far 
use only one color at a time. / or i? bands are 
typically chosen for shear measurement because 
of the increased galaxy number density, advanced 
detector technology, and small PSF at these wave- 
lengths. Surveys like C0MB017 (Brown et al. 
2003), and VIRMOS/Descartes (van Waerbeke et 
al. 2002) are leading a trend to use additional 
multicolor photometry to provide photometric 
redshifts of the source galaxy population. The 
SNAP surveys will be simultaneously observed in 
9 bands: 6 optical colors spanning roughly B ^ I, 
plus J, H' and K (the near IR filters are twice as 
large and receive double the total exposure times 
given in §2.1; see paper I). We have not simu- 
lated this multicolor data, btit it will inevitably 
raise the S/N of shear estimation for every source 
galaxy. At a minimum, image coaddition or si- 
multaneous fits to shapes in several colors will 
increase the effective exposure time. Something 
more ambitious, like shifting to the rest-frame R 
or the rotating disc dis-alignment suggested by 
Blain (2002), might even reduce systematic mea- 
surement biases. Further work is needed in cosmic 
shear methodology to investigate the optimal use 
of multicolor data. However, it can already be 
said that our current monochromatic approach 
will yield a conservative estimate of the lensing 
sensitivity expected from future analyses. 

3. Weak Lensing sensitivity 

In this section, we determine the accuracy with 
which it is possible to recover the input shear from 
the noisy image simulations. The formalism of 
shapelets can be used to form an accurate shear 
measurement (Refregier & Bacon, 2003). How- 
ever, since the images themselves were created us- 
ing shapelets, we choose here to be conservative 
and use a slightly older but independent method 
developed by Rhodes, Refregier & Groth (2000; 
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RRG). 

3.1. Advantages of space 

We first discuss tlie advantages specific to weak 

lensing mcasiircmcnt that arc provided by obser- 
vations from space. Tiie figure of merit for any 
lensing survey needs to include more than the 
etendue, a product of the survey area and the flux 
gathering power of a telescope (Tyson et al. 2002, 
Kaiser et al. 2002). It must also account for the 
finite PSF size, the size-magnitude distribution of 
background galaxies, and systematics {e.g. due to 
the atmosphere or telescope optics). Shear sensi- 
tivity is raised for a spacecraft over a ground-based 
telescope for the additional reasons listed below. 

• More objects have measurable shapes. Al- 
though not as much sky area will be sur- 
veyed as by proposed ground based sur- 
veys such as MEGACAM (Boulade et al. 
2000), VISTA (http://www.vista.ac.uk), 
or LSST (http://www.lsst.org), the num- 
ber density of resolved objects is an order 
of magnitude higher from space (compare 
figure 5 with Bacon et al. 2001). Such an 
increase in S/N per unit area will enable 
the mapping of projected dark matter maps 
with adequate resolution for a direct com- 
parison to redshift surveys (§5.1) and the 
generation of a mass-selected cluster cata- 
log [e.g. Weinberg & Kamionkowski 2002, 
Hoekstra 2002). Quantitative study of high- 
sigma mass fluctuations is one of the most 
promising methods to break degeneracies in 
cosmological parameter estimation, particu- 
larly constraining flm {e-fl- van Waerbeke & 
Mellier 1997; Cooray, Hu & Miralda-Escude 
2000; Munshi & Jain 2001; Schneider 2002). 
Furthermore, studying well-resolved groups 
and clusters individually, rather than statis- 
tically, will lead to a better understanding of 
astrophysical phenomena such as biasing or 
the mass-tc!rnpc!rature relation (Weinberg ct 
al. 2002; Huterer & White; 2003 Smith et al. 
2003). 

• The shape of individual galaxies are more 
precisely measured. The SNAP PSF is small 
(0.13" FWHM assuming 4/xm CCD diffu- 
sion). It is more isotropic and, importantly, 
more stable than even the HST PSF (see 



paper I). This enables shape measurement 
to be more reliable, or possible at all, for 
small, distant galaxies. The stable photome- 
try from the 3-day orbit may even permit the 
use of weak lensing magnification as well as 
shear information (see e.g. Jain 2002; 2003). 
Whether directly measured or inferred from 
shear, this in turn is useful to correct for 
the effect of lensing on the distance moduli 
to the SNAP supernovas (Dalai et al. 2003; 
Perlmutter et al. 2002). 

• Galaxy redshifts are known accurately and 
to a greater depth. SNAP's stable 9-band 
optical and NIR imaging is ideal to produce 
exquisite photometric redshifts for almost all 
galaxies at 2; < 3 detected at 5ij in the I- 
band (see §4). This should be compared to 
the ~38% completeness of photo-2s possi- 
ble from the ground in the COMBO-17 data 
with a similar cut and a median redshift of 
R c± 2A (Brown et al. 2002). This allows 
a good estimation of the redshift distribu- 
tion of source galaxies, the uncertainty in 
which is a major contribution to the error 
budget in current lensing surveys. Projected 
2D power spectra and maps can be drawn in 
several redshift slices, using redshift tomog- 
raphy. More ambitiously, cluster catalogs 
and dark matter maps can be constructed 
directly in 3D (§5.2), enabling the 3D corre- 
lation of mass and light and the tracing of 
the growth of mass structures. 

• Galaxies are farther away. Distant objects, 
too faint and too small to be seen from the 
ground, are measurable from space. The 
evolution of structures can thus be traced 
from earlier epochs, giving a better handle 
on cosmological parameters (see paper III). 
Furthermore, recent numerical simulations 
(Jing 2002; Hui & Zhang 2002) suggest that 
intrinsic galaxy alignments impact lensing 
surveys to a greater depth in redshift than 
previously assumed. If this is confirmed, in- 
trinsic alignments will mimic and bias cos- 
mic shear signal in all but the deepest sur- 
veys, where the galaxies are farther apart 
in real space. Using 3D positions of galax- 
ies from SNAP photo-zs, it will be possi- 
ble to isolate close galaxy pairs and to mea- 



7 



sure their alignments, or to optimally clown- 
weight close pairs thus reducing their impact 
(King & Schneider 2002; Heavens & Hey- 
mans 2002). 

3.2. Shear measurement method 

The advantages of space-based data described 
above will provide limited gains without an 
equally precise and robust shape measurement 
method. The now-standard weak lensing method 
for ground-based data was introduced by Kaiser, 
Squires & Broadhurst (1995; KSB). KSB forms 
shear estimators from quadrupole and octopole 
moments of an object's flux. Modern techniques 
are being developed to incorporate higher order 
shape moments or Baycsian statistics to raise 
the sensitivity to shear. These methods include 
shapelets (Refregier & Bacon, 2003) and oth- 
ers by Bernstein & Jarvis (2001); Bridle et al. 
(2003); Kaiser (2000). However, since the simula- 
tions themselves were created using shapelets, we 
choose here to be conservative and use the inde- 
pendent method developed by Rhodes, Refregier 
& Groth (2000; RRG). This is related to KSB, 
but optimized for use with space-based data. It 
has already been used extensively on HST images 
(Rhodes et al. 2001; Refregier et al. 2002) and is 
therefore appropriate for our current purposes. 

Following KSB, RRG measures a galaxy's 
two-component ellipticity £i from the Gaus- 
sian weighted quadrupole moments of its surface 
brightness I (9), 



where 



J a 



^ _ {Jii — J22, 2J12} 
Jii + J22 

J <fe Oj 0j w{e) i{9) 



(1) 



(2) 



and w{9) is a Gaussian of width adjusted to match 
the galaxy size. The unweighted PSF moments arc 
measured from a (simulated) starfield and RRG 
corrects the galaxy ellipticities to first order for 
PSF smearing. Occasional unphysical ellipticities, 
\e\ > 2, are excluded, along with galaxies fainter 
than AB 26.5 (for the SNAP wide survey) or AB 
29.1 (for the SNAP deep survey) and with sizes 



R: 



Note that R is an rms size measure rather than a 
FWHM, and that this procedure does indeed se- 
lect only resolved objects. The locations of these 
cuts have been chosen to yield reasonably stable 
results; the effect of moving the size cut is dis- 
cussed further in section §3.4. 

RRG finally provides the shear susceptibility 
conversion factor, G, to generate unbiased shear 
estimators 7^ for an ensemble of objects, given by 



7i = 



M 
G 



(4) 



where G depends upon the fourth order moments 

Jijki of a ga-la-xy population, defined similarly to 
equation [2]. In our simulated SNAP images, G is 
of order 1.6. 
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Fig. 4. — The applied shear 7i„ in the shapelet 
simulated images vs its recovery ^out using an inde- 
pendent measurement method (Rhodes, Refregier & 
Groth 2000). The image used is one 7.5 arcmin^ re- 
alization of the SNAP deep survey shown in figure 1. 
The recovery is linear, but the slope of the fitted line 
(dotted) is flatter than that expected (dashed line). 

This shear measurement method and the simu- 
lations are tested in figure 4. An artificial shear is 
applied uniformly upon all objects in a 7.5 arcmin^ 
simulated image, in the 72 = and 71 = di- 
rections, before convolution with the SNAP PSF. 
Using RRG, we correct for the PSF smearing and 
recover the input shear. As can be seen in figure 
4, the recovery is linear, but the slope (see dotted 
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line) is underestimated (dashed line). This incon- 
sistency probably has two origins: inaccuracy of 
the image simulations and instabilities in the shear 
measurement method. The latter may be removed 
with future techniques. For the purposes of this 
paper, we follow the procedure adopted by Bacon 
et al. (2001), where a similar bias was observed 
in the KSB method. We apply a linear correction 
factor to the measured shears and to their errors. 
This factor is (0.79 ± 0.03)"^ at the depth of the 
HDF, and (0.87 ± 0.04)"! for the SNAP wide sur- 
vey. 

Even after this correction, there remains a small 
difference in the rms scatter of galaxy ellipticities 
between the simulations and real Groth strip data 
(Rhodes, Refregier & Groth 2001). As shown in 
Massey et al. (2003), this discrepancy is not de- 
tected with the standard shape measures of SEx- 
TRACTOR (Bertin & Arnouts 1996), however RRG 
proves to be a more sensitive test. Perhaps be- 
cause of the precise properties of the background 
noise, or perhaps because the wings of simulated 
objects are truncated beyond the SExtractor 
isophotal cutoffs, is observed by RRG to be 
lower in the simulated images by another factor of 
~ 0.8. Work is in progress to establish the origin 
of this effect. For the purposes of this paper, we 
simply increase the error bars by this amount. 

3.3. Shear sensitivity of SNAP 

Now that the image simulation and analysis 
pipeline is in place, we can measure SNAP's sen- 
sitivity to shear. Trade-off studies are under way 
for several alternative telescope designs, including 
the level of CCD charge diffusion, the pixel size, 
the effect of DRlZZLEing, and the coefficient of 
thermal expansion in the secondary struts, which 
may be the main cause of temporal variation in the 
PSF (see paper I). Here we present the results of a 
study which uses the baseline design specifications 
and time-averaged PSF of the SNAP satellite. In 
this study the PSF used is the residual between 
the HST and SNAP PSFs (see §2.3), in order to 
keep the size distribution of galaxies realistic for 
SNAP images. 

The top panel of figure 5 shows the surface num- 
ber density rig of galaxies in a survey of a given 
exposure time texp on SNAP. The exposure times 
reflect a ^ 5 x overall improvement in instrument 
throughput and detector efflciency over WFPC2 
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Fig. 5. — Shear sensitivity as a function of SNAP 

exposure time tcxp- Top panel: the surface number 
density of all galaxies (rig, tot) detected by SExtrac- 
tor and of the subset (rig) of these useable for weak 
lensing, i.e. having survived further cuts in size and 
ellipticity by RRG (see text). Second panel: the me- 
dian / band magnitude, Jm, in the two subsets of the 
galcixy catalog, which has been interpreted as median 
redshift, Zm, using equation [8]. Third panel: the rms 
error cr-^ = (lll^)^''^ per galaxy for measuring the shear 
7, after PSF correction and shear calibration. Bottom 
panel: the rms error a— for measuring the mean shear 
7 in 1 arcmin^ bins. The dot-dashed line shows an es- 
timate of the expected rms shear in a ACDM universe. 
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on HST (Lampton ct al. 2002). The dashed hne 
shows the number density rig, tot of all the galaxies 
detected by S Extractor, after a S/N cut which 
is equivalent to /< 29.1 at the depth of the HDF. 
As discussed in §3.2, galaxies which are too faint, 
too small, or too elliptical are excluded from weak 
shear catalogs. The solid line shows the number 
density Hg of galaxies which are useable for weak 
lensing following the magnitude, size and elliptic- 
ity cuts. The error bars reflect the uncertainty in 
measuring number counts at low S/N and an es- 
timated sample variance between the HDF-N and 
HDF-S. 

An important cut in the weak lensing analysis 
is the size cut, which reduces the detected galaxy 
sample by about 30% at the depth of the HDF. 
This fraction is a strong function of PSF size, and 
is thus much larger for ground based imaging. As 
can be inferred from the top panel of figure 5, the 
SNAP wide survey {ug ~ 100 galaxies arcmin"^) 
will thus provide a dramatic improvement over 
current ground-based surveys {Ug ~ 25 galaxies 
arcmin"^ are used by most groups; see e.g. Bacon 
et al. 2002). The effect of moving the size cut is 
discussed further in section §3.4. 

The second panel of figure 5 shows the median 
magnitude, Im, of the galaxy catalog before and 
after cuts in size and ellipticity by the weak lens- 
ing analysis software. This has been converted to 
a median redshift, z^, using equation [8]. For the 
purposes of this plot, we assume that this relation- 
ship is still valid even after the size cut. 

The third panel of figure 5 shows the rms er- 
ror (T-y = (|7p)^^^ per galaxy for measuring the 
shear, after the PSF correction and shear calibra- 
tion. The slightly increasing error at longer texp 
reflects the decreasing size of fainter galaxies, and 
correspondingly less resolved information content 
available about their shapes. To map the shear, 
the noise can be reduced by binning the galaxies 
into cells. The rms noise of the shear 7 averaged 
in a cell of solid angle A — 1 arcmin'^ is given by 



and is plotted in the bottom panel of figure 5. The 
wide and deep SNAP surveys will thus afford a la 
sensitivity for the shear of ~ 3.0% and better than 
2.2% on this scale, respectively. As a comparison, 
the rms shear expected from lensing on this scale 



in a ACDM model is approximately 3% (assum- 
ing 0^=0.3, nA=0.7, (78=0.9, r=0.21). This sig- 
nal increases with survey depth because the total 
lensing along a line of sight is cumulative. The 
wide SNAP survey will thus be ideal to map the 
mass fluctuations on scales of 1 arcmin^, with an 
average S/N of unity in each cell. The recovery of 
simulated mass maps will be discussed in §5. 

Note that the shear sensitivities presented here 
arc conservative estimates, particularly for the 
deep SNAP survey. The image simulations ex- 
tend so far only to the depth of the HDF. Future 
shear measurement methodology will also be more 
accurate and stable on any individual, resolved 
galaxy than the RRG method used in this pa- 
per. Higher order shape statistics {e.g. shapelets) 
will be used, as will simultaneous measurements 
in multiple colors and pre-selection of early-type 
galaxy morphologies. 

3.4. Effect of size cut and pixel scale 

Small, faint and highly elliptical objects are ex- 
cluded from the final galaxy catalog in the RRG 
shear measurement method. Of all these cuts, it is 
the size cut that excludes the most objects. In an 
image at the depth of the HDF, about 30% of de- 
tected galaxies are smaller than our adopted size 
cut at R = 1.7 pixels. The exact position of this 
cut has been determined empirically to produce 
stable results, from experience with both HST 
data and our simulated images. The quantitative 
effects of moving the size cut are demonstrated in 
figure 6. 

If the cut is moved to a larger size, fewer objects 
are allowed into the final galaxy catalog, and the 
shear field is sampled in fewer locations. Conse- 
quently, both dark matter maps and cosmic shear 
statistics become more noisy. If smaller galaxies 
are included in the catalog, the shear field is indeed 
better sampled, but the shape measurement error 
is worse on these galaxies. The bottom panel of 
figure 6 shows that moving the size cut to smaller 
objects has no net change in the precision of shear 
recovery: adding noisy shear estimators to the cat- 
alog neither improves nor worsens the measure- 
ment. A size cut at R = 1.7 pixels is optimal at 
the depth of the HDF and in the observing condi- 
tions modelled by our image simulations. To sim- 
plify comparisons of galaxy number density, the 
same cut has been applied to data at the depth of 
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Fig. 6. — Shear sensitivity as a function of size cut 
R in the RRG sliapo measurement metiiod for the 
wide SNAP survey (sohd line) and at the depth of 
the HDF (dashed line). The vertical dotted line shows 
the fiducial value adopted elsewhere in our analysis. 
Top panel: the surface number density of galaxies use- 
able for weak lensing. Middle panel: the rms error 
a~f = {I'fl^y^^ per galaxy for measuring the shear 7, 
after PSF correction and shear calibration. Bottom 
panel: the rms error for measuring the mean shear 
7 in 1 arcmin^ bins. 

the SNAP wide survey. A different cut could have 

been adopted, producing fewer galaxies but each 
with more accurate shear estimators: the crucial 
figure ax/ would not change. (This is especially 
true in the SNAP wide survey because of the rel- 
ative dearth of small galaxies). 

As described in section 2.4, we have assumed 
that an effective image resolution of 0.04" can 
be recovered for SNAP data by taking multi- 
ple, dithered exposures, and either stacking them 
with the DRIZZLE algorithm or by fitting each 
galaxy's shape simultaneously in them all. The in- 
crease in image resolution from these techniques is 
vital for cosmic shear measurements. The number 
density of useable galaxies increases dramatically. 



and the measurement of their shapes is improved. 
Were it not possible to apply DRIZZLE or to re- 
cover this resolution, the large pixel scale currently 
proposed for SNAP would seriously impair shear 
measurement. A size cut at i? = 0.12" (= 3 pixels 
in figure 6) would roughly halve the number den- 
sity of useable sources and correspondingly reduce 
the sensitivity to gravitational lensing. 

4. Photometric Redshift Accuracy 

Gravitational lensing is achromatic, so shear 
measurement may be performed in any color. 
As discussed in §2.4, current techniques measure 
galaxy shapes in only one band at a time (usu- 
ally R 01 I are chosen for their steeper slope of 
number counts). However, gravitational lensing 
is also a purely geometrical effect, and measure- 
ments are aided greatly by accurately knowing 
the distances to sources. The latest surveys, and 
future high-precision measurements will therefore 
require multiple colors for photometric redshift 
(photo- 2;) estimation. Reliable photo-2;s will not 
only remove current errors due to uncertainty in 
the redshift distribution of background sources, 
but even make possible an entirely 3D mass recon- 
struction, as demonstrated in §5.2, Taylor (2003a), 
Hu & Keeton (2002), Bacon & Taylor (2002) and 
Jain & Taylor (2003). 

snap's thermally stable, 3-day long orbit is 

specifically designed for excellent photometry on 
supernovse. Combining all 9 broad-band filters (6 
optical, 3 NIR) will also provide an unprecedented 
level of photo-z accuracy, for all morphological 
types of galaxies over a large range of redshifts. In 
this section, we simulate SNAP photometric data 
in order to determine this precision. 

We have used the hyperz code (Bolzonella, Mi- 
ralles & Pello 2000) to generate the observed mag- 
nitudes of a realistic catalog of galaxies following 
Lilly et al. (1995), 

^(/) ^ 10°-^^x^ , (6) 

where / is the /-band magnitude. The galaxies 
were assigned a distribution of Spectral Energy 
Distribution (SED) types similar to that in real 
data and containing ellipticals, spirals and star- 
burst galaxies. Redshifts were assigned at random, 
and independently of spectral type, according to 



11 



Koo ct al. (1996) as verified by the DEEP collab- 
oration (1999), 



SNAP wide 



SNAP deep 



dN 
dz 



{z) ~ z'e 



where 



0.722 + 0.149(1-22.0) 



(7) 



(8) 



(Lanzetta, Yahil & Fernandez- Soto 1996). SNAP 
colors were then inferred by integrating the SED 
across filter profiles, adding an amount of noise 
corresponding to the exposure time and instru- 
ment throughput. 

Hyperz was then used again, to estimate red- 
shifts for the simulated catalog as if it were real 
data. Unlike the image simulations in §2, this ap- 
proach can already be taken to the depth of both 
the wide and the deep SNAP surveys by extrapo- 
lating functional forms for the luminosity and red- 
shift distributions (eqs. [6] & [7]). Magnitude cuts 
were applied at AB 26.5 (wide) or AB 29.1 (deep) 
in R. Similar magnitude cuts were made in each 
filter, chosen at the lOa detection level of an ex- 
ponential disc galaxy with FWHM=0.12" (Kim et 
al. 2002). Past experience with lensing data (see 
§3.2, Bacon et al. 2000) confirms that this is rea- 
sonable S/N limit. Note however that the size 
and cUipticity cuts implemented for the simulated 
images in §3.2 were not included at this stage. 

Figure 7 shows the precision of photometric red- 
shifts in both the wide and deep SNAP surveys. 
All galaxy morphological types are included in this 
analysis. Clearly demonstrated is the need for the 
near IR HgCdTe detectors, a component of the 
satellite where a spacecraft has a clear advantage 
over the ground. Figure 8 shows the accuracy of 
the photo function of source (photomet- 

ric) redshift. Here, Azpijoto(-z) is the rms of the 
core Gaussian in a double-Gaussian fit to horizon- 
tal slices through the distributions in figure 7. 

To estimate the accuracy of 3D mass recon- 
structions (§5.2), we now concentrate on objects 
closer than z ~ I. According to equations [6] and 
[7], these make up ~ 38% of all galaxies detected 
in R for the wide SNAP survey, and ~ 35% for 
the deep. For the lensing analysis (§3.2), we have 
to reject some fraction of galaxies because they 
were too small & not resolved. Here, we assume 
that the same percentage of rejection applies to 




Fig. 7. — Recovery of redshifts of a realistic popula- 
tion of galaxies using hyperz with the SNAP filter set. 

Top-left panel: the wide survey, using all 9 colors. Top- 
right panel: the deep survey using all 9 colors. Bottom 
row: the same, but with only the 6 optical colors, as 
if the near IR HgCdTe data were not available. 



the z < 1 sub-sample of galaxies. This will yield a 
conservative estimate of the number of objects re- 
maining in the real SNAP survey because objects 
closer than z=l are likely to have a larger median 
size than the entire sample. Removing this frac- 
tion from the number density of galaxies shown in 
figure 5 leaves 40 ± 5 useful galaxies per arcmin^ 
in the SNAP wide survey to 2; = 1 , and more than 
90arcmin~^ in the SNAP deep survey. For these 
galaxies only, Azdeep =0.034 and Azu,i(ie =0.38 
using all nine SNAP colors. 
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Fig. 8. — Azphoto, the rms scatter on photometric er- 
ror estimation, as a function of increasing source red- 
shift, z. Top panel: results for the SNAP deep survey. 
Bottom panel: results for the SNAP wide survey. In 
both cases, the solid line shows photometric rcdshift 
errors using observations in all 9 SNAP optical and 
near IR colors. The dashed lines show errors if the 
near IR HgCdTe data were unavailable. 

5. Dark Matter Mapping 

In this section, we describe the prospects of a 
space-based weak lensing survey for mapping the 
2D and 3D distribution of dark matter. Because 
of the high number density of background galaxies 
resolved from space, this is one application where 
a mission like SNAP will fare particularly better 
than surveys from the ground. 

5.1. 2D Maps 

To simulate observational data, we begin with 
shear maps created by raytracing through N-body 
simulations from Jain, Seljak & White (2000). 
We then add noise to these idealized data, cor- 
responding to the predicted levels for SNAP or 
observing conditions at the currently most suc- 
cessful ground-based facilities. In each case, wc 
then attempt to recover the input projected mass 
distribution by inverting the map of the shear into 
a map of the convergence k. Convergence is pro- 
portional to the projected mass along the line of 



sight, by a factor depending on the geometrical 
distances between the observer, source and lensed 
galaxies (see e.g. Bartelmann & Schneider, 2001). 

Figures 9 and 10 show how the projected mass 
can be mapped from space and from the ground. 
The color scale shows the convergence k. The top 
panel of figure 9 shows a (noise-free) simulated 
convergence map from the ray tracing simulations 
of Jain, Seljak & White (2000) for an SCDM 
model. Underneath it is a version smoothed by 
a Gaussian kernel with an rms of 1' for compari- 
son to the simulated recovery from observational 
data in figure 10. 

Figure 10 shows similarly smoothed mass maps 
that would be possible using (from top to bot- 
tom) a ground-based survey, the SNAP wide sur- 
vey and the SNAP deep survey. These were pro- 
duced by adding to k, before smoothing, Gaus- 
sian random noise to each 1 arcmin^ cell with an 
rms of (T— given by equation [5] . Overlaid contours 
show mass concentrations detected at the 3a, 4cr 
and 5(T levels. For ground based observations, we 
set (7^ = 0.39, and used Ug = 25 arcmin"^, as is 
available for ground-based surveys {e.g. Bacon et 
al. 2002). For the SNAP wide and deep survey, 
the surface density of useable galaxies was taken 
to be Ug = 105 arcmin^^ and 259 arcmin^^, and 
the rms shear noise per galaxy was taken to be 
(j-y =0.31 and 0.36, respectively, as derived from 
figure 5. The galaxies are assumed to all have a 
redshift of z=l, which is a good approximation 
as long as the median redshift is approximately 
that value. As noted above, the surface density 
and median redshift will actually be higher for the 
SNAP deep survey, because only exposure times 
corresponding to that of the HDF were simulated. 

From the ground,only for the strongest features 
(i.e. the most massive clusters) can a 3(7 detec- 
tion be obtained. From space, the very high den- 
sity of resolved background galaxies allows the re- 
covery of uniquely detailed maps, including some 
of the filamentary structure and individual mass 
overdensities down to the scale of galaxy groups 
and clusters. Thus, SNAP offers the potential of 
mapping dark matter over very large fields of view, 
with a precision well beyond that achievable with 
ground-based facilities. 

The masses and locations of individual clusters 
can be extracted from such maps using, for ex- 
ample, the Map statistic (Schneider 1996, 2002), 
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Fig. 9. — Left column. 2-dimensional convergence 
maps projected along a lino of sight. The convergence 
K is proportional to the total matter density along the 
line of sight, and can be deduced from the shear field. 
Top panel: simulated (noise-free) convergence map de- 
rived by ray-tracing through an SCDM N-body sim- 
ulations of large-scale structure from Jain, Soljak & 
White (2000). The region shown is 30'x30' and the 
sources are assumed to lie at z = 1. Bottom panel: 
same map but smoothed using a Gaussian kernel with 
a FWHM of 1'. 

Fig. 10. — Right column. Reconstructions of the con- 
vergence map in figure 9 which may be feasible from 
weak lensing surveys on the ground and from space. 
Overlaid contours show 3a, 4a and 5a detection lim- 
its. Top panel: convergence k with noise added corre- 
sponding to ground-based observations {i.e. Ug — 25 
arcmin"^ and a-y = 0.39; Bacon et al. 2002). Mid- 
dle panel: convergence map with the expected noise 
properties of the wide SNAP survey {i.e. rig = 105 
arcmin"^ and cr-y = 0.31). Bottom panel: with the 
expected noise level of the deep SNAP survey {i.e. 
Ug = 259 arcmin"^ and a^ = 0.36). 
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which has been apphed successfully to find mass 
peaks in several surveys {e.g. Hoekstra et al. 
2002a; Erben et al. 2000) and our own work; or 
the inversion method of Kaiser & Squires (1993; 
KS) which was used by Miyazaki et al. (2002). 
Marshall, Hobson & Slosar (2003) demonstrate 
the effectiveness of maximum entropy techniques 
to identify structures in KS lensing maps, using 
criteria set by Bayesian evidence. White et al. 
(2002) argue that using any detection method, a 
complete mass-selected cluster catalogue from 2D 
lensing data would require a high rate of false- 
positive detections, since the prior probability 
is for them to be anywhere throughout a given 
survey. This has been avoided in practice by 
secondary cross-checks of the lensing data with 
spectroscopy, deep x-ray temperature or SZ ob- 
servations. Indeed, two previously unknown clus- 
ters have already been found in weak lensing maps 
and spectroscopically confirmed by Wittman et al. 
(2001, 2003). However, this confusion does make 
it harder to resolve the debate on the possible ex- 
istence of baryon-poor "dark clusters" {e.g. Dahle 
et al. 2003). These are a speculative population of 
clusters which would be physically different to and 
absent from the catalogues of optically or x-ray 
selected clusters. Remaining dark- lens candidates 
(Erben et al. 2000; Umetsu & Futamase 2000; Mi- 
ralles et al. 2002) have currently been eliminated 
as chance alignments of background galaxies (or 
possibly associations with nearby ordinary clus- 
ters Gray et al. 2001; Erben et al. 2003). If others 
could be found in high S/N weak lensing maps, 
they would present a challenge to current models 
of structure formation, and need to be accounted 
for in estimates of flm', but they would be unique 
laboratories to decipher the nature of dark matter. 

5.2. 3D Maps 

The growth of mass structures can be followed 

in a rudimentary way via photometric redshifts, 
by making 2D mass maps or power spectra with 
source galaxies in different redshift slices (see e.g. 
Bartelmann & Schneider 2000 §4) . This technique 
is useful for a global statistical analysis of a survey 
in order to constrain cosmological parameters. It 
is used as such in paper III, to predict possible con- 
straints with SNAP. Tomographic measurements 
of shear have also led to estimates of mass and 
radial position of clusters (Wittman et al. 2001, 



2003). After this analysis, spectroscopic redshifts 
were needed to constrain the mass further by fix- 
ing the precise radial position of clusters. 

An alternative approach, in which one naturally 
reclaims the radial mass information as well as the 
transverse density, has been developed by Taylor 
(2003a) and Hu & Keeton (2002). In this method, 
the shear pattern on an image is treated as a fully 
3D field, by including from the outset the redshift 
of galaxy shear estimators as well as their 2D po- 
sition on the sky. Taylor (2003a) shows that there 
is a simple inversion that relates this 3D distortion 
field to the underlying 3D gravitational potential. 

Using this technique, we now demonstrate the 
capabilities of SNAP for reconstructing the 3D 
mass distribution and locating clusters. We ap- 
ply the simulations of Bacon & Taylor (2003) to 
the telescope and survey parameters deduced in 
paper I, §2 and §3 then attempt to recover the 
gravitational potential of two M = IO^^Mq NEW 
clusters at redshifts of 0.25 and 0.4 and separated 
by 0.2 degrees on the sky (see Bacon & Taylor, fig- 
ures 4 & 5). Note that this is specifically a search 
for clusters, which induce a significant shear signal 
at one location, rather that integrating the impact 
of many small objects and filamentary structures 
in a statistical basis over an entire shear field. Our 
relatively simple input model is therefore appro- 
priate for our current purposes: it is a common 
occurrence that, for a line of sight with large shear, 
a single cluster along a line of sight is responsible 
for the signal. 

Eirst, we calculate the corresponding lensing 
potential for this field (using the prescription of 
Bacon & Taylor, equations 9-12; c.f. Kaiser & 
Squires 1993). As the lensing potential field 
is an integral of the shear field, we are able to 
reconstruct (j) with more accuracy than the grav- 
itational potential which is a function of the 
second derivative of (p; nevertheless (p itself con- 
tains valuable 3D information. This is discussed 
in full in Bacon & Taylor (2003). We have used the 
expected number density of useable galaxies closer 
than z = 1 to be a conservative ng{z < 1) = 90 
per arcmin^, combining results from figure 5 and 
§4. In this nearby regime we can approximate 
^ oc z^. We have taken into account the shot 

dz 

noise arising from intrinsic galaxy ellipticities, us- 
ing an error on shear estimators for galaxies of 
a-y =0.36. We have also included an error on our 
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photometric redshifts of Azphoto = 0.034 through- 
out < z < 1, from §4. 

Figure 11 shows the reconstruction of the lens- 
ing potential out to z = 1 available with the SNAP 
deep survey. The units of the lensing potential 
here are radians^, having chosen the differential in 
K = 0.5d^<f> to be taken in units of radians. In this 
simulation, we see that the lower redshift cluster 
is very pronounced in the lensing potential, with 
S/N per pixel of 5.4 at z—1. The lensing potential 
due to the higher redshift cluster is also clearly 
visible, with S/N per pixel of 3.0. 

Figure 12 shows a reconstruction of the 3D 
gravitational potential, using Taylor's inversion 
and Wiener filtering (Bacon & Taylor equations 8, 
40). Even with a simulated mass of only lO^^M©, 
the lower redshift cluster is very pronounced, and 
the higher redshift cluster is also detectable at 
the 4.4(7 level. Extra noise peaks in figure 12 
demonstrate that the extremely low end of cluster 
catalogues will be subject to high false detection 
rates. However, this reconstruction affords mea- 
surement of masses of matter concentrations to an 
accuracy of AM ~ 1.1 x lO^^M© at ^ ~ 0.25 or 
AM ~ 4.8 X 1O^^M0 at z ~ 0.5 via fitting (c./. 
Wittman et al. 2001, 2003). We can also estimate 
radial position of mass concentrations from the 
simulated lensing data with accuracy Az ~ 0.05 
for clusters of mass M = IO^^Mq at z < 0.5 (c./. 
Bacon & Taylor 2003). The mass concentrations 
are observed at z = 0.25 and 0.4 as expected with 
peak S/N of 2.8 and 3.3 respectively (N.B. this 
is S/N per pixel; the overall detection significance 
of the cluster is as quoted above). Of course, the 
sensitivity of this technique drops for clusters at 
greater distances, as their induced lensing poten- 
tial grows less within the observed redshift window 
(i.e. z < 1). In an alternative regime of interest, 
mass fluctuations J ~ 1 are measurable on degree 
scales (c./. Hu & Keeton 2002). 

Equivalent simulations can be carried out for 
ground-based experiments, providing more limited 
prospects. The key difference that makes a space- 
based experiment superior over a ground-based ex- 
periment in this regard is the reduced error on 
shear estimates for galaxies, particularly for galax- 
ies at z > 0.5, due to improved resolution and 
small PSF. From the ground, studies of the 3-D (j) 
field are restricted to measuring the mass of a clus- 
ter along the fine of sight at the AM ~ 2 x lO^^M© 
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Fig. 11. — Top panel: reconstructed lensing potential 
from a finite number of z < 1 galaxies with realistic 
ellipticities; Ug = 90 arcmin , (^^^ — 0.36, as expected 
for the SNAP deep survey. The a;-axis represents an- 
gle in degrees; the j/-axis represents redshift. The two 
crosses mark the positions of the input M = 10^'' M0 
NFW clusters. Bottom panel: difference between in- 
put and recovered lensing potential fields. 
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Fig. 12. — Reconstructed gravitational potential us- 
ing galaxy properties of the SNAP deep survey, as in 
figure 11. The a-axis represents angle in degrees; the 
y-axis represents redshift. The two input clusters are 
clearly seen as the two darkest regions near the center 
of the image, marked with a cross. 

level a,t z = 0.25, with a 1.3a measurement of a 
mass of a.t z = 0.5 along the line of sight 

of a foreground z = 0.25 cluster. Reconstruction 
of $ in 3D is possible on 5' scales only out to a 
redshift of ^ ~ 0.5 (see Bacon & Taylor 2003). 
Application of the full 3D inversion technique to 
real ground-based data is currently being carried 
out, and even measurements of one cluster behind 
another cluster are possible (see review by Taylor 
2003a). 

snap's ability to measure the 3D gravitational 
potential in this fashion is of great importance. 
One can determine the mass and density profile of 
several matter concentrations along a line of sight, 
avoiding the ambiguity of surface-density lensing 
or projection effects, and obtaining accurate mea- 
surements of the mass of matter clumps in 3D. 
One can directly compare the visible matter dis- 
tribution with the underlying mass distribution 
to obtain important information regarding bias- 
ing and galaxy formation as a function of redshift. 
One can also examine the number of objects ex- 
ceeding a certain mass threshold as a function of 
redshift (see e.g. Viana & Liddle 1996), or recon- 
struct the 3D power spectrum directly (see Taylor 



2003a) in order to obtain constraints on cosmologi- 
cal parameters or to test the gravitational instabil- 
ity paradigm which is thought to govern structure 
formation. 

6. Conclusions 

We have shown how a space-based wide field 
imager like SNAP is ideally suited for studies of 
weak gravitational lensing. The aspects of this 
satellite's design relevant for weak lensing, and the 
baseline survey strategy, were presented in paper 
I. A shapelet-based method for creating simulated 
space-based images (Massey et al. 2003, Refregier 
2003) has been used to predict SNAP's sensitivity 
to shear, taking explicitly into account its instru- 
mental throughput, limitations and sensitivity. 

In this paper, we have considered the baseline 
SNAP design for our predictions. As explained 
in paper 1, this design is almost optimal because 
many requirements to find supernovae are the same 
as those to measure weak lensing (stable imaging, 
small PSF, excellent multicolor photometry). 

The increased image resolution available from 
space makes possible the construction of high res- 
olution projected dark matter maps with an rms 
shear sensitivity of ~2.5% in every 1' cell for the 
300 square degrees wide SNAP survey and better 
than 1.8% for the SNAP deep survey (c./. expected 
mean signal in a ACDM universe is approximately 
3%). Since lensing is sensitive to mass regard- 
less of its nature and state, these maps will be 
unique tools for both astrophysics and cosmolog- 
ical parameter estimation. Statistical properties 
of the dark matter distribution will be precisely 
measured at several cosmological epochs and con- 
straints on Om, (Tg and w are discussed in paper 
III. 

SNAP'S simultaneous 9-band observations also 
open up new opportunities for 3-dimensional map- 
ping via photometric redshift estimation (Taylor 
2003a, Hu fc Keeton 2002, Bacon & Taylor 2002). 
snap's photometry allows an excellent resolution 
of Az = 0.034 in redshift. Here we have shown 
that SNAP will measure mass concentrations in 
full 3D with a la sensitivity of approximately 
1 X 1O"M0 at ^ ~ 0.25 and ~ 5 x IO^^Mq at 
2: ~ 0.5. In this fashion it will be possible to di- 
rectly trace the non-linear growth of mass struc- 
tures, testing with high precision the gravitational 
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instability theory. 

Space-based wide-field imaging can be com- 
bined with weak gravitational lensing to produce 
2D and 3D mass-selected cluster catalogs down to 
the scale of galaxy groups. Mass and light in the 
local universe can be mapped out with exquisite 
prccdsion, thus offering exciting prospects for both 
astrophysics and cosmology. 
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